o 

(N 



Lie algebra automorphisms as Lie point symmetries and the 
solution space for Bianchi Type /, 77, IV, V vacuum geometries 

Petros A. Terzis* and T. Christodoulakis^ 

University of Athens, Physics Department 

Nuclear & Particle Physics Section 

Panepistimioupolis, Ilisia GR 157-71, Athens, Hellas 



Abstract 



£^-i Lie group symmetry analysis for systems of coupled, nonlinear ordinary dif- 

ferential equations is performed in order to obtain the entire solution space to 
Einstein's field equations for vacuum Bianchi spacetime geometries. The symme- 
tries used are the automorphisms of the Lie algebra of the corresponding three- 
dimensional isometry group acting on the hyper-surfaces of simultaneity for each 
Bianchi Type, as well as the scaling and the time reparametrization symmetry. 
\^q \ The method is applied to Bianchi Types I, II, IV and V. The result is the acqui- 

sition, in each case, of the entire solution space of either Lorenzian of Euclidean 
signature. This includes all the known solutions for each Type and the general 
solution of Type IV (in terms of sixth Painleve transcendent Pvi)- 
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1 Introduction 

The exploitation of the group of automorphisms in order to obtain a unified treatment of 
spatially homogeneous Bianchi spacetime geometries has a rather long history as it goes 
back to the early 1960's [lj. In 1979, Harvey [2] found the automorphisms of all three- 
dimensional Lie algebras, while the corresponding results for the four- dimensional Lie 
algebras have been presented in [3J. In Jantzen's tangent space approach the automor- 
phism matrices are considered as the means for achieving a convenient parametrization 
of a full scale factor matrix in terms of a desired, diagonal matrix jU [5j [6]. Siklos used 
these time-dependent automorphisms as a tool for the proper choice of variables aiming 
at a simplification of the ensuing equations [7], while Samuel and Ashtekar were the first 
to look upon automorphisms from a space viewpoint [5]. The notion of Time- Dependent 
Automorphism Inducing Diffeomorphisms ,i.e., coordinate transformations mixing space 
and time in the new spatial coordinates and inducing automorphic motions on the scale- 
factor matrix, the lapse, and the shift vector has been developed in [9]. The use of these 
covariances enables one to set the shift vector to zero without destroying manifest spa- 
tial homogeneity. At this stage one can use the rigid automorphisms, i.e., the remaining 
gauge symmetry, as Lie point symmetries of Einstein's field equations in order to reduce 
the order of these equations and ultimately completely integrate them [10] . 

The methodology applied in the present work is the Lie group symmetry analysis 
of differential equations. For a detailed account of this approach see, e.g., [UJ, [12], 
[T5] . Just to put it simply, a symmetry of a system of differential equations is a trans- 
formation mapping of any solution of the system to another solution. In other words, 
the symmetry group of the system is the group transforming solutions of the system to 
other solutions. In fact, it is the largest local group of transformations acting on the 
system variables. Such groups are Lie groups depending on continuous parameters and 
consisting of either point transformations, acting on the systems space of independent 
and dependent variables or, more generally, contact transformations that act also on all 
first derivatives of the dependent variables. The theory can be also generalized to the 
so called higher-order (Lie-Backlund) symmetries, by including in the transformation of 
the independent variables derivatives of the dependent. In this work we use all these 
three kinds. In general, the aforementioned transformations are nonlinear and the main 
benefit of the Lie group symmetry analysis comes from the replacement of the nonlinear 
symmetry conditions of the system investigated by "linear" conditions associated with 
the infinitesimal generators of the symmetries. To this purpose, it is necessary first to 
determine a class of general admissible variable transformations and then to search for 
special members of this class under which the system of differential equations remains 
invariant. It is rather evident, that the degree of generality of the admissible transforma- 
tions is proportional to the number of the existing symmetries. One should further stress, 
that one of the main ingredients of the method lies on the notion of "prolongation" of 
a group action on the space of derivatives of the systems dependent variables up to any 
finite order. Thus one is able to deal with differential equations of any order. In closing 
this short introduction, we must point out that the knowledge of a symmetry group of a 



higher order differential equation has much the same consequences as the knowledge of 
a symmetry group of the corresponding system of first order differential equations. In 
fact, we will make use of the following two important theorems (for a proof, see [12]) : 

Theorem 1.1 Let 

— = F v (x,y), u = l,...,q (1.1) 

be a first order system of q differential equations and suppose that G is a one-parameter 
symmetry group of the system. Then, there exists a change of variables (t,u) = ip(x,y) 
under which the system can be written as 

du u 

— = H v (t,u\...,v?- 1 ), v = l,...,q (1.2) 

so that the original system is reduced to a new system of q — 1 differential equations for 
u 1 ,!! 2 , ,u q ~ 1 together with the integral 



u q 



(t) = J H q (t, u\t),..., u q -\t)) dt + c. (1.3) 



Theorem 1.2 If the system of q first order differential equations considered in Theo- 
rem 1.1 admits an r-parameter solvable symmetry group, then the solution of the system 
can be found by quadrature from the solution of a reduced system of q — r first order 
differential equations. If the original system is invariant under a q-parameter solvable 
symmetry group, then its general solution can be found by quadrature alone. 

The paper is organized as follows: in Section 2, we give our method. In Sections 3-6, 
the application of the method to Bianchi Types I, II, V and IV is presented. Finally, 
some discussion and concluding remarks are given in Section 7. 

2 The Method 

It is known, that the line element for spatially homogeneous spacetime geometries with 
a simply transitive action of the corresponding isometry group |14j . [15], assumes the 
form 

with the 1-forms cr a defined from 



ds 2 = (N a N a - N 2 ) dt 2 + 2N a a?dx i dt + - fa i 3 ay^dx i dx j (2. 1) 



d<r a = C^y A *< & o% - a*, = 2C^a7a?. (2.2) 

(Small Latin letters denote world space indices while small Greek letters count the dif- 
ferent basis one-forms; both types of indices range over the values of 1, 2, 3). 
Then the Einstein field equations are (see, e.g. [16], [9]): 

E = K afi K afi - K 2 - R = (2.3) 



E a = KZC; t - K'C^ = (2.4) 

E aP = K aP + N (2K r a K r/3 - KK afS ) + 

+ 2N» (K av C v Pp + K 0V CQ - iVR Q/3 = 
where 

K a/3 = ~ (7^ + 2 lav C^ p N» + 2 7/3i ,C^iV) (2.6) 

is the extrinsic curvature of the three-dimensional hypersurface and 

R a/3 = C^C^y^ypyf^ + 2CfyC*, + 2C^ A7|a/7 ^+ 



2C^C^ 7aA7 - + 2C^C^fixlT 



(2.7) 



is the Ricci tensor. 

In [9], particular spacetime coordinate transformations have been found, which reveal 
as symmetries of (J2.3J) . (12. 4p . and (12. 5p the following induced transformations of the 
dependent variables N, N a , 7a/3 : 

N = N, N a = K{N p + lpa P°) ) ~ lpv = kl^ vlafi (2.8) 

where the matrix A and the triplet P a must satisfy: 

K C M = C ^ A "^ (2.9a) 

2P»C^ = A" (2.9b) 

These transformations were first presented in jl], see also the discussion on p. 3586 of 

0- 

For all Bianchi Types, this system of equations admits solutions that contain three 

arbitrary functions of time plus several constants depending on the automorphism group 
of each Type. The three functions of time are distributed among A and P a (which also 
contains derivatives of these functions) . So one can use this freedom either to simplify 
the form of the scale factor matrix or to set the shift vector N a to zero. The second 
action can always be taken, since, for every Bianchi Type, all three functions appear in 

pa 

In this work we adopt the latter point of view. Having used the freedom stemming 
from the three arbitrary functions in order to set the shift vector to zero, there is still a 
remaining "gauge" freedom consisting of constant A^ (automorphism group matrices of 
the Lie group defined by the structure constants C| ' . Indeed, the system ( 12 .9p accepts 
the solution A» = const., P a = 0. The latter are also known in the literature as rigid 
symmetries [T7] . 

The generators of the corresponding motions ^ pv = A p A^ 7a/ g induced in the space 
of the dependent variables spanned by 7^'s (the lapse is given in terms of 7a/ g, %,£ by 
algebraically solving the quadratic constraint equation), are: 

Xj = \ p Ialp p^- (2.10) 

O^a/3 



with A^ satisfying 

tfp C m = K/3 ^7 + ^% C 1p- ( 2 - ! !) 

These generators define a Lie algebra and each one of them induces, through its 
integral curves, a transformation on the configuration space spanned by the 7 Q /3's [18J. 
If a generator is brought to its normal form (i.e., ^-), then the Einstein field equations, 
written in terms of the new dependent variables, will not explicitly involve 2$. They thus 
become a first order system in the function ij [llj. If the aforesaid Lie algebra is abelian, 
then all generators can be brought to their normal form simultaneously. If the Lie algebra 
is non-abelian, then we can diagonalize in one step those generators corresponding to 
any eventual abelian subgroup. The rest of the generators (not brought in their normal 
form) continue to define a symmetry of the reduced system of the field equations if the 
Lie algebra of the X^'s is solvable [12] • One can thus repeat the previous step by 
choosing one of these remaining generators and bring it to its normal form. This choice 
will of course depend upon the simplifications brought to the system at the previous 
level. Finally, if the Lie algebra does not contain any abelian subgroup, one can always 
choose one of the generators, bring it to its normal form, reduce the system of equations, 
and search for its symmetries (if there are any). Lastly, two further symmetries of 
(12. 3p . (12 ,4p . and (I2.5p are also present and can be used in conjunction with the constant 
automorphisms: The time reparametrization t — > t + a, owing to the non-appearance of 
time in these equations (the system being autonomous), and the scaling by a constant 
7«/3 - > A7 a/ 3 (homothety) as can be straightforwardly verified. Hence, in every Bianchi 
Type there are, added to the Xm generators, also the following generators: 

r, = I (2.12) 

d d d d d d 

07ll O712 0713 d7 22 ^723 Ol33 

These generators commute among themselves, as well as with the Xm's: 

[Xj,Y a ]=0 {/ = l,2,3,4|a = l,2} (2.14) 

3 Bianchi Type I 

In the Bianchi Type I model, the structure constants are 

C a o = for every value of a, (3, 7 (3.1) 

By using this relation and the definition of the 1-forms cr a (|2.2p . we get 

cr 1 = dx, cr 2 = dy, a 3 = dz (3.2) 

with the corresponding Killing fields 

fi = d x , 6 = d y , £3 = d z (3.3) 



From equations (I2.9p we have that the triplet P a (t) is arbitrary while the automorphisms 
A? are constant, i.e. they are elements of the general linear group GL 3 (R). Hence, the 
symmetry group of the field equations is a ten-parameter group consisting of the nine 
variables of the GL3(R) group plus the symmetry produced by the generator (I2.12p . 
while the generator (I2.13P is already included in the GL 3 (M.) group. Using the triplet 
P a (t), we can set the shift vector N a to zero by applying (12.81) . Furthermore, in order 
to simplify the field equations, we choose the time gauge iV 2 = 7, 7 = det(j a p). In this 
case, eq. (j2.4p is satisfied identically, while equations (I2.3p . (J2.5P take the form 

1^ l Vfi laB % V - 1^ T la? %u = (3.4) 

lad ~ 1 PT jpa IrB = (3.5) 

The last equation can be integrated if multiplied by 7°°": 

i a<T %p-i aa i PT % a <jr(] = 0^ 

7 a<T 7 a /3 + 7 aCT 7 pr 7p a 7^ = 0^ 
1™ 7^+7^7^ = 0^ 

7 - 7q/3 = ^, 0*, = const. (3.6) 

Multiplying the latter by j ff p we obtain the linear system 

7/3p = $° 3 lap or, in matrix form, 7 = ^7 (3.7) 

Thus, the quadratic constraint (13.41) can be written as a relation between the constant 
elements of the matrix 1? 

tfV^a - CO 2 = ° ( 3 - 8 ) 

The general solution of the linear system (13.71) is of the form 

j a3 = exp ($ T t) p a c pB , exp ($ T t) p a c pB = exp (§ T t) P ^ c pa , c a8 = const. (3.9) 

where the symmetry of 7^ has been taken into consideration. 

In order to find the analytic form of this solution, it is necessary to calculate the 
matrix exponential of the 3x3 matrix $. To this purpose, we first note that 7^ 
transforms, under GL 3 (R), as 

7 ^=A^A^7^^7 = A T 7A, A E GL 3 (R) (3.10) 

so we have for the linear system (13.71) 

A T ^A = ^ T A T 7 A^^ = ^ 7; ^=A" T ^ T A T (3.11) 

Now, A is an arbitrary matrix, hence it can be used to simplify the matrix d. The degree 
of simplification depends on the eigenvalues of the matrix d which may be (i) three real, 
distinct eigenvalues, (ii) one real and two complex conjugate eigenvalues, (iii) three real 
eigenvalues two of which are repeated and, (iv) three real, repeated eigenvalues. If the 
matrix d is the zero matrix, then ^ a Q = const, and the metric thus obtained describes 
Minkowski flat spacetime . 



3.1 Three real distinct eigenvalues 

In this case, we can choose the matrix A such as to diagonalize the matrix $, i.e. 

& = diag(p 1 ,p 2 ,P3) => exptfi = diag(e Plt ,e P2 \e P3t ). (3.12) 

As a result, the solution obtained from (13. 9 j) is 

7 = diag( Cl e plt } c 2 e P2 \ c 3 e P3t ) (3.13) 

and the constants (pi,p 2 ,P3) are related through (13.81) . i.e. 

Pl+Pl+Pl = (Pi +P2 +P3) 2 =^PlP2 +PlP3+P2P3 = (3.14) 

Since the constants (^1,^2,^3) are different, none of them vanishes. Thus, dividing by, 
let's say, pi, we can eliminate it, so 

Pl + PJL + PlPJL = ^ a + p + a p = , a = P Aft = ^ (3.15) 

Pi Pi Pi Pi Pi Pi 

The coefficients of the diagonal elements of the matrix 7 can, by a proper coordinate 
transformation of (x, y, z), become the same and equal to p± 2 . By putting t = p^ 1 r, we 
get the following final form of the metric: 

ds 2 = - e {1+a+ ^ T dr 2 + e T dx 2 + e aT dy 2 + e? T dz 2 (3.16) 

where the constants (a, ft) satisfy (13.15J) . The multiplicative constant p{ 2 does not 
appear in the metric due to the homothety 

H = 2 (a + 1) d T + a 2 x d x + y d y + (a + l) 2 z d z . (3.17) 



The above metric was first given, although in a different form, by Kasner |19j . The 
metric is particularly interesting for the values (a, ft) = (1, —1/2) or (a, ft) = ( — 1/2,1): 
In addition to the three Killing fields f)3.3p . there is a fourth of the form 

U = yd x -xd y . (3.18) 

The non permitted pair of values (a, ft) = (0, 0) leads the metric (I3.16P to the standard 
Minkowski form, so we can include these values into the domain of the constants (a, ft). 

3.2 One real and two complex conjugate eigenvalues 

In this case, the matrix i? cannot be diagonalized, but we can choose an apprpriate A 
such as to bring i? into its "rational normal form", i.e., 

'pi \ 

$= ( 1 ] (3.19) 

-{p 2 2 +pD 2p 2 J 



with eigenvalues (pi , p 2 + i p 3 , P2 — i Pz) ■ Thus, we have for the exponential matrix exp ($ t) 



/e pi4 \ 

g*P2 gP2* 

(-p 2 sin(p 3 t) + p 3 cos(p 3 t)) sinp 3 t 

P3 P3 

gP2* g*P2 

ip\ + p\) sin(p 3 1) (p 2 sin(p 3 i) + p 3 cos(p 3 *)) 



(3.20) 



and the constants (pi,P2,P3) are related through (|3.8p : 

2pip 2 +P2+P3 = (3.21) 

Since the constants pi, p 3 cannot vanish, we can divide the previous relation by pfp 3 
obtaining 

2A + /3A 2 + /3 = 0^/3 = --^, /3 = ^,A = ^ (3.22) 

A 2 + 1 px p 3 

Utilizing the fact that the matrix j a p must be symmetric as well as the freedom of using 
the automorphisms, we get, after the transformation t — (3pit, the following metric: 

ds 2 = e (2X + ^)r dT 2 + e r/(] dx 2 _ e Xr sinTd 2 + 

X X 2 ( 3 - 23 ) 

+ 2 e T cos r dy dz + e T sin rdz 

which possesses a homothety produced by the field 

H = -4Xd T -4X 2 xd x + (y (-A 2 + 1) + 2 \z) d y - (z (A 2 - 1) + 2 \y) d z (3.24) 

This metric was first given, although produced in a different way, by Harrison [201 • Also 
in this case there are special values of the constant A for which we have a fourth Killing 
field: 

/3 

A = ^-^£ 4 = 6<9 r + 2 V3xd x - (\/3 y + 3 z) d y + (3 y - V3 z) d z 

/3 
A = -^- =► & = 6<9 T - 2 y/3xd x + (V3y -3z)d y + (3y + V3z) d z , 

while there is no homothety. Finally, it is worth noting that in metric f !3.23|) the hyper- 
surface t = const, is spacelike. 

3.3 Three real (two repeated) eigenvalues 

In this case, the matrix $ can be brought in its Jordan normal form by a proper choice 
of the matrix A, i.e. 




gpit 











gP2t 


e P2 H 








e P2t 



= p 2 1 =^exp(tf£) = e P2t e P2 *t (3.25) 



where the constants (pi,P2) are related through (13. 8ft : 

P2(P2 + 2pi) = (3.26) 

Here too, the constant p\ is not zero, so dividing the previous relation by p\ yields 

A (A + 2) = ^ A = V A = -2, — = A (3.27) 

P\ 

With the help of the automorphisms and the transformation r = p\t we get the metric 

d s 2 = e (2X+1)r dr 2 + e T dx 2 + t e Xr dy 2 + 2e Xr dy d z (3.28) 

This metric is also a member of the Harrison class [20J and admits a homothety produced 
by the field 

H = 2 d T + 2 x A d x + (A + 1) y d y + (z (A + 1) - y) d z (3.29) 

In this case too, the hypersurface t = const, is spacelike. Furthermore, for the value 
A = the metric (13.281) describes a pp-wave, since the Killing field u = £3 = d z has zero 
covariant derivative and zero measure: 

\ = 0^u a u a = 0A u a .p = (3.30) 

3.4 Three real repeated eigenvalues 

In this case, the Jordan normal form of the matrix $ is 

. /M0\ _ (It \t 2 \ 

0= p 1 ^ exp(tft) = e p M 1 t , (3.31) 

\0 p) \0 1 / 

so (13. 8p gives p = 0. Again, exploiting the fact that the matrix 7^ is symmetric as well 
the automorphism matrices, we get the metric 

ds 2 = dt 2 + 2t 2 dx 2 + dy 2 - Atdxdy + Adxdz (3.32) 

which has the following properties: 

• There is a homothety, with corresponding field 

H = td t + yd y + 2zd z (3.33) 

• It is a pp-wave metric since, for the Killing field u = £3 = d z , we have 

u a u a = A u a .p = Q (3.34) 



• The hypersurface t = const, is spacelike. 

At this point, the question arises whether the two metrics (13.281) with A = and (13.321) are 
actually the same since they both are pp-wave metrics, they have three Killing fields one 
of which is timelike, and they both admit a homothety. The classical way to answer this 
question is by turning to the scalar invariants that can be constructed by the Riemann 
tensor and its covariant derivatives. This approach is not applicable to our case, since 
the fact that the metrics are pp-wave metrics implies the vanishing of all these scalar 
invariants. An alternative approach is to find a tensor that is identically zero for one 
of the metrics and non zero for the other. This would mean that there is no coordinate 
transformation relating the two metrics. If we consider the tensor 

^■a^Se — R K a j3 -Rft7<5e;A, (3.35) 

we observe that it vanishes identically for metric (I3.32p but not for metric (I3.28P with 
A = 0. Thus, the two metrics describe different geometries. 

4 Bianchi Type II 

In this case, the structure constants are 
r* 1 — —r x — 2 

23 ~~ ° 32 — 9 

C a g = for all the other values of a, (3, 7 (4.1) 

Utilizing these values and the definition (12. 2p of the 1-forms cr a , we get: 



'z 1 0\ 

1 \ ^ a- 1 = zdx + dy, cr 2 = dz, a 3 = dx (4.2) 

1 0/ 

with the corresponding Killing fields 

fi = dx, 6 = -xd y + d z , £ 3 = d y (4.3) 

From (12. 9p we have for the automorphisms A a g and the triplett P a (t) 

y i+s2 -s 3 s 4 x(t) y(ty 

A a =\ e si s 3 j (1.1a) 

s 4 e S2 



P a (t) = [ P(t), e ?~ g3 " , g4 T e ^ ) (4.4b) 

1 w Vi +S2 -S3S4 e Sl+S2 -S3S4 / 

10 




where P(t), x(t) and y(t) are the three aforementioned arbitrary functions of time which 
appear in the vector P a and can be thus used to set the shift vector N a to zero. The 
remaining symmetry of the field equations is described by the constant matrix 



3 Sl+S2 _ 



M c 








S3 S 4 S5 S$ 

e Sl s 3 

s 4 e 



This matrix defines a transformation of the variables j a 0, 7^ 
be written analytically as 



(4.5) 



M^ a M v ^^ v , which can 



711 

7l2 
713 



D Sl+S2 



S1+S2 



D S!+S 2 



- s 3 s 4 ; 711 

- s 3 s 4 ) (e Sl 712 + s 4 713 + s 5 711) 

- s 3 s 4 ) (e S2 713 + s 3 712 + s 6 711) 



„2si 



722 = e Sl 722 + 2 e Sl s 4 723 + 2 e Sl s 5 712 + s 4 733 + 2 s 4 S5713 + s 5 711 

723 = S3 (e si 722 + s 4 723 + s 5 712) + e S2 (e Sl 723 + s 4 733 + s 5 713) + 

s 6 (e Sl 712 + s 4 713 + s 5 711) 



(4.6) 



k 733 = e S2 733 + s 3 722 + 2 s 3 s 6 712 + s 6 711 + 2 e S2 (s 3 723 + s 6 713) 

These equations define a transformation group G r of dimension r = dim(Aut(II)) = 6. 
The six group generators are determined through the definition 



X, 



ds A 



d 



s=0 



•97, 



afi 



(4.7) 



where A = {1,2,3,4,5,6}. Application of this relation to (14. 6p results in the following 
explicit form of the six generators: 

X 1 = 2 711 <9 711 + 2 712 d 112 + 713 <9 713 + 2 7 22 <9 722 + 723 9 723 

X 2 = 2 711 <9 7ll + 712 <9 7l2 + 2 713 <9 7l3 + 723 9 723 + 2 733 <9 733 

A 3 = 712 <9 7l3 + 722 d j23 + 2 723 <9 733 

X A = 713 <9 712 + 2 723 <9 722 + 733 d 123 

X 5 = 711 ^712 + 2 712 <9 722 + 713 <9 723 

A 6 = 711 <9 7l3 + 712 <9 723 + 2 713 <9 733 (4.8) 

The algebra g r associated with the group G r has the following table of non-vanishing 
commutators: 



[*i,*s] 

[X 2 ,X 3 ] 

[A 3 ,X 4 ] 



x 3 
-x 3 

X\ — X2 



[Xi,X 4 ] 

[x 2 ,x 4 ] 

[X 3 ,X 5 ] 



-x A 

x A 

-x 6 



[Xi,X 6 ] 


= x 6 


[x 2 ,x 5 ] 


= x 5 


[X^Xs] 


= -x 5 



(4.9) 



As it is implied by the previous commutators (14.91) . the group is non abelian, hence we 
cannot simultaneously diagonalize all the generators. However, the generators X 5 ,X 6 



11 



Starting 


«6 — "5 ^ 

from the 


u. 

equations 


E 2 = 


0,^3 


= 0, we 


g 


variables ii 2 . 


u 3 : 
















E 2 '- 


= 0^ 


U5U2 - 


-u 4 ii 3 = 


= 






E 3 = 


= 0^ 


u 6 u 2 - 


■ W 5 U 3 = 


0. 



commute with each other, as well as with the generator Y 2 , so we can bring them into 
their normal form by the following transformation of the 7 a/ g: 

711 = e Ul(4) , 712 = e u ^u 2 (t), 7l3 = e Ul « u 3 (t) 

I22 = e u ^(u 2 2 (t)+u,(t)) (4.10) 

723 = e Ul(t) (u 2 (t) u 3 (t) + u 5 (t)), 733 = e Ul(t) (^(t) + « 6 (i)) 

► The solution space in the gauge N a = 

Before we proceed to the solution of the Einstein equations, we have to find the values 
allowed for the functions Ui, % — 1, . . . , 6. The determinant of the matrix 7 a/ g, is 

detfrcp] = e 3 " 1 ( Ui u 6 -uf) , (4.11) 

0, we get the following system for the 

(4.12a) 
(4.12b) 

This system admits non zero solutions only if u^Uq — u\ = 0, a condition forbidden by 
our restrictions. Thus, we have 

u 2 (t) = k 2 , u 3 (t) = k 3 (4.13) 

Now, by using the automorphism matrix 

(I -k 2 -k 3 \ 
M a p = 1 , (4.14) 

\0 I ) 

the matrix 7^ is brought to block diagonal form, and we have the following theorem: 
The general solution of Bianchi Type II corresponds to a block diagonal form in the basis 
of the 1 -forms cr a . 

As we have now proved that the matrix 7 a/ g is in block diagonal form with no loss of 
generality, we can repeat the former procedure for the calculation of the generators, i.e., 
we can start from the matrix 

/7n \ 

7a/3 = I 722 723 (4.15) 

V 723 733/ 

with the corresponding automorphism matrix 

si+s 2 _ g 3 g 4 

M a =- j e Sl s 3 | . (4.16) 

s 4 

12 




The generators thus obtained are 

X x = 2 7u <9 711 + 2 722 <9 722 + 723 <9 7 23 
X 2 = 2 7n <9 7ll + 723 <9 723 + 2 733 <9 733 

^3 = 722 <9 723 + 2 723 <9 733 

X 4 = 2 723 <9 722 + 733 <9 723 . (4.17) 

The generators Y 2 , Xi + X 2 , X 3 commute with each other, so they can be brought to 
their normal form by applying the transformation 

7a/? = | e Ul+2u6 u 4 e Ul+2u6 u 4 u 5 |. (4.18) 

e Ul+2u6 u 4 u 5 & 

In this parameterization, the only variable of which the second derivative appears in the 
field equations is w 4 , since the generators have been transformed to 

Y 2 = 7 —,X 3 = —,X 1 + X 2 = — (4.19) 

OU\ ou 5 ou e 

To further simplify the present situation, we shall also use the generator Yj_ = jj?. The 
variable u 4 , that does not appear in the generators Y 2 , X\ + X 2 , X 3 , can be utilized in 
order to parameterize time t, i.e., we can apply the transformation 

t ->. U4 (t) = s, u x (t) -¥ ui{s), u 5 (t) -> u 5 (s), u 6 {t) -> u 6 {s) (4.20) 

Of course, this transformation can be applied only in the case where u±(t) is not a 
constant function, hence, we must distinguish between two cases: either the function 
u±{t) is constant or not. 

4.1 The case u^t) = const. 

Let u 4 (t) = k 4 . Then, the linear constraints (12. 4 p are identically satisfied, while the 
quadratic constraint (12. 3p defines the lapse function N 2 (t) as 

N 2 = k A e Ul (3 ill -hii 2 5 + 16 u x u 6 + 20 iif) (4.21) 

Substituting the latter into the dynamic equations (12.51) . we can solve them for the 
variables U\, tic,, Uq. Starting from equation En = 0, we notice that the coefficient of ii 5 
is proportional to the quantity 

u 5 (mi + 4m 6 ). (4.22) 

Therefore, in order to solve for u^, we must be sure that this coefficient is not zero. The 
case Mi + 4it 6 = leads to a zero lapse, so only the case u${t) = k 5 is left. 
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Assuming that u^(t) 7^ k 5 , we solve equation En = for u 5 and, substituting in 
equation E 2 2 = 0, we can solve the latter for U\. The coefficient of U\ is the quantity 
k^iiQ which cannot vanish since, for k± = 0, we have det(^ a p) = while, for u§(t) = fc 6 , 
we obtain a zero lapse. Now, substituting in E23 = we get 

fc 4 ii 5 (iii + 3 uq) = =>• U\ = k\ — 3 u 6 . (4.23) 

However, this leads to 1x5 = ^5, which is impossible in this "branch" of the solution. 

If we start with w 5 = fc 5 , then the coefficient of iii in equation £^n = is proportional 
to the quantity 

u 6 (mi + 3m 6 ) (4.24) 

The case u G = k 6 leads to a zero lapse, so the case remaining to be considered, before 
we solve equation En = for the variable u\ , is 

u 6 (t) = h- ^. (4.25) 

With this value for the function u 6 (t) the rest of the equations are satisfied identically, 
but the lapse function takes now the value 

N i = A e Mt) kiil j, (4.26) 

9 

i.e., we have ended up with a metric of Euclidean signature. Using the automorphism 
matrix 



M 



/Vhe~ 2k(i \ 

^k A k\ + 1 e~ k& 



q __k±k§__ — kg V&4 g~ fc 6 1 

\ y/k 4 kj + l y/k 4 fef+1 / 



(4.27) 



the matrix 7 Q/ 3 can be brought in diagonal form. By choosing a gauge in which ui(t) = 
3 lnt, we get the metric 

d s 2 = t dt 2 + - (a 1 ) 2 + t (a 2 ) 2 + t (a 3 ) 2 (4.28) 

In addition to the three Killing fields (14. 3p . this metric has a fourth one: 

& = -2zd x + (z 2 -x 2 )d y + 2xd z (4.29) 

and it admits a homothety produced by the field 

H = td t + xd x + 2yd y + z d z (4.30) 

Since we have examined the case in which the coefficient of iii in En = vanishes, we 
can solve this equation for this function, so 

u t = — (3 ii\ + 25 ill ug + (68 ii 2 6 + u 6 ) u x + 60 iif) , (4.31) 

Uq 
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hence the rest of the equations are also satisfied. In order to simplify the latter equation 
we choose a gauge in which u§(t) = 3t, thus obtaining the autonomous first order 
equation for iii 



ill = u\ + 25 u\ + 204 tij + 540 

This equation can be integrated with the help of the transformation 

1 



t = y(f)> «i = ^ «i 
which brings the equation in the following form: 



y'(0 



(4.32) 



(4.33) 



1 



(£ + 10)(£ + 9)(£ + 6) 



|/(0 = fcj + \ Hi + 10) - ~ ln(£ + 9) + 1 ln(£ + 6) 



(4.34) 



Using the automorphism matrix 



/ £tj/3 _6fci-fc 2 /2 



M 



V2 



\ 



V 



o 





_27 „-3fci 

2 C 



|/xfc 5 e- 3fel - fc2 / 2 



2 e -3fci-fc 2 /2 



(4.35) 



and redefining the constant k± and time £ as 



8\/3 



k 4 = r-Ll^g-fca^ 



^ = -6 



5e 2,/, -l 



(4.36) 



27 r " * " 3e 2 ^-l 

we obtain the metric 

ds 2 =fi 2 (-e 2 ^ cosh V> dt 2 + sech^ (cr 1 ) 2 + e^ cosh^ (<r 2 ) 2 + e^ cosh^(o- 3 ) 2 ) (4.37) 

In addition to the three Killing fields (I4.3p . this metric has a fourth one: 

^4 = -2 29 x + (2 2 -x 2 )^ + 2x9 2 (4.38) 

and it does not admit a homothety, so the constant /x cannot be absorbed. 

4.2 The case u^t) ^ const. 



Let us assume that u±(t) = t. Then the linear constraints ( 12. 4 p are satisfied identically, 
while the quadratic constraint (I2.3P defines the following lapse function N 2 (t): 



N 2 = e Ul (3tul - t 2 ii 2 5 + 2 (8tu 6 + 1)^ + 2 (WtiiQ + 3) u 6 ) 



(4.39) 
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Substituting this in equation E u = 0, we find that the coefficient of u.5 is proportional 
to the quantity M5 (u\ + 4m 6 ) which cannot vanish since it leads to a zero lapse. Solving 
equation En = for u 5 and then substituting in equation E 2 2 = 0, we find that the 
coefficient of uq is proportional to the quantity tiii + 2 which cannot vanish because in 
such a case we obtain again a zero lapse. Thus, we can solve for u§ and substitute in the 
rest of the equations. Equation £"23 = includes the function iti the coefficient of which 
is proportional to the quantity 

6tu 6 + 2tux + l, (4.40) 

the vanishing of which leads, as we shall see, to a metric of Euclidean signature. Assuming 
that the quantity (I4.40p is not zero, we get the following system of differential equations: 

fi 5 = -^(t 3 u\ + 1) (4.41a) 

u s = — (-6tul + t 2 ul{-4tu 6 + 3) -4w 1 (8£w 6 + 1) - 2w 6 (20£w 6 + 7)) (4.41b) 

Ul = - (ptul - 5t 2 ii 2 5 + Qu 6 (Wtu 6 + 3) + (48tu 6 - 2t 3 ii 2 5 + 5) Wi) (4.41c) 
We shall first consider the case of a vanishing quantity (14.40J) : 

6 1 m 6 + 2 1 iii + 1 = => u 6 = k 6 In t u x (4.42) 

6 3 

Solving equation Eu = for u 5 and substituting in E22 = 0, we obtain 

m 5 = => u 5 = k 5 (4.43) 

that brings equation E 33 = into the form 

3£ 2 m 1 + £ 2 m? + 7£mi + 4 = 0. (4.44) 

The latter is a Riccati equation for w\(t) = iii(t) 

^ = -\ wl --L Wl -jr^ (445) 

which is easily solvable since a partial solution of it is already known: w\ = — |. With 
the transformation 

o 1 

(4.46) 







Wt(t) = - 


2 

-7 + 


1 
h(t) 


equation 


(I4.45I) 


takes the linear form 







h= -h + -=>h(t) = - Int + kit. (4.47) 

1(3 



By integrating, we finally have for Ui(t) 

Ux = k 2 -2 lnt + 3 ln(lnt + 3fci). 
Redefining the constants k 2 , k§ 

k 2 = ln/z 2 — 3 ki, 
using the automorphism matrix 



kf, 
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(3 lnK 2 -ln/? + 3£;i), 



(4.48) 



(4.49) 



/. /Z p -3fci/4 



M 



\ 



e 








. 4 /£fi e 3fci/8 





VK e -9fei/8 



(4.50) 



/ 



. 4 /— e -9 ^/ 8 

and putting £ = e^ _3fcl , we obtain the metric of Euclidean signature 

ds> = /i 2 (e-t £de + -^ {a 1 ? + Z (<x 2 ) 2 + e -« £ (<x 3 ) 2 ) • (4.51) 

This metric does not admit a homothety, so the constant jjl cannot be absorbed. 

Now, returning to the system f !4.4ip and multiplying the first equation by W5, we have 



u 5 u 5 

y 



2i4 (fi4 + K 



t 

Ay 



dul 
~dt 



Av 2 
^ul + l) 



(t 3 y+l), y = ii 



5- 



(4.52) 



i.e. a Ricatti equation for y{t). A partial solution of this equation is y\ = —473, hence 
we can reduce it to a linear differential equation by applying the transformation 



1 1 
+ 



At 3 h 



and get 



h = - h + At 2 =>- h = --r^t 2 + At' 6 => u b = k 52 ± \lk 51 . 



(4.53) 



(4.54) 



t &5i " ' V ""* t 

Consequently, the last two equations of the system (14.411) become 

tit = (iii\ Ax\u e ), u 6 = (u 1 \A 2 \ii 6 ), (4.55) 

where use has been made of the notation (ui\ = (liiiiif) and \ui) = (««!' with the 3x3 
matrices Ai, A 2 given by the relations 

5 li 



/ 



-4-i 



60 



At 2 (k 5l t-1) t 

10Afcit-ll 

48 



2k^t 2 -2t 



9 



(4.56) 



0/ 
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Ao 



14 fat -13 2Q \ 



V 



;t 2 (A; 51 t-l) 2k 51 t 2 -2t 

-16 

t 
-3 / 



(4.57) 



The latter system is of polynomial form in the variables ii\, Uq, hence we can simplify it 
based on the transformation 



1 . 1 

4 (fc 51 1 - 1) 4t 



Thus, we obtain the system 

4 2/1 2/2 + hi 



Vi 



Ay/t(k 51 t-lY 



y-2 



i.e., 



2/2 = 2>i =^ 2/2 = 2/i + fr V fc. 



l/i 



42/12/2 + fcgi 

4 ^ (fri * - 1; 



42/ 2 + 4fc 2 v / fri2/i + fri 



4v/t(fc 5 it-l) 
The previous Riccati equation is simplified by applying the transformation 



t H- -: — cosh (2r), 2/1 (£) h-> vfrij/ifr), fc 2 i-> -cosh/i, 
K51 



thus arriving at the equation 



1 



2/i + 42/i — 4 cosh /i 2/1 + 1 = 



2/1 = - [cosh/z — sinh/x tanh (2 (fci — r) sinh/x)] . 

Substituting 2/1 (r) in (I4.58p . we obtain 

Mi —C\ —At cosh/i + 3 ln(cosh(2 [k\ — r) sinh/i)) + In — 

Uq =c 2 + r cosh u — ln(cosh(2 (fr — r) sinh a)) — -In — — 

4 smh4r 



(4.58a) 
(4.58b) 

(4.59) 



(4.60) 



(4.61) 



(4.62) 

(4.63a) 
(4.63b) 



In order to get the final form of the metric, we redefine the constants \x = 2 a, k 5 i = 
exp(ci — Aki cosh2a)/2 k 2 , we choose the gauge r = k\ — £csch2a/2 and we use the 
automorphism matrices 



M 



' — ne 2 1 2c2 


V ° 








\ 

Jf2 e~ c 2+ fc i- fc i cosh 2 o- 



K — J — C2 — fci+fci cosh2cr k — J — C2+k\+ki cosh2cr 
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for /c 52 = a/&5i and 



with 



' d\ d 2 — d%di 
.} / = | d x d 3 

<i 4 d 2 . 



e% +V2Kh 2 e 2k i™^° __a_ ( ,__, | _, j< „^^ 



d 1 = — - — — g-^-C2+fcl-fci coshia (4.64) 



— 9"— C2— fci+fci cosh 2<t 



d 2 =- i7 = e~ 2 - C2 - fel+fel coshi(T (4.65) 

-rf _ /o ^ U „2fci cosh 2(7 

rf 3 = " V2ft: ^ e e -f -*-*i-*i ^2. (4 66) 



4 = ^= e -f-«+fei+fci cosh2<x ^ 467 ^ 

for A)52 7^ 1//C51. Thus, we arrive at the form 

ds 2 = K 2 fe^ coth2a cosh^+sech^f+e^ ^ cosh£ (<r 2 ) 2 +e« tanhff cosh£ (<r 3 ) 2 ) 

(4.68) 
This solution was first produced by Taub [21]. It does not admit a homothety, hence 
the constant k cannot be absorbed. It is noteworthy, that in the limiting case a — > +oo, 
this metric reduces to metric ( I4.37P which possesses four Killing fields. 

5 Bianchi Type V 

In this case, the structure constants are: 

r l — —H l — r 2 — —r 2 — - 

u 13 ~~ u 31 ~~ u 23 ~~ ° 32 ~~ r, (5 1) 

C a g = for all the other values of a, (3, 7 

Using these values in the definition (14 .7p of the 1-forms er a we have 

/0 e~ x \ 
a a { = e~ x ^ a 1 = e~ x dz, a 2 = e~ x dy, <r 3 = dx (5.2) 

\1 0) 

with corresponding Killing fields 

£1 = d z , 6 = d v , & = d x + yd y + z d z (5.3) 

The time- dependent automorphisms are: 

f Pl p{t) P2 p{t) x (ty 

A a = I p 3 P(t) p 4 P(t) y(t) I , Pl p 4 - p 2 p 3 = 1 (5.4) 

1 
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and 



P a (t) = x{t) In 



xjt) 
Pit) 



,y(t) 



Pit)) 



(lnP(t))' 



(5.5) 



where P{t), x(t) and y(t) are arbitrary functions of time that can be utilized in order to 
set the shift vector to zero. The remaining symmetry of the field equations is described 
by the constant matrix: 

fe Sl s 2 s 3 \ 
M = s 4 e S5 s 6 (5.6) 

\0 l) 

where we have changed the parameterization of the constants so that the choice Sj = 
corresponds to the identity matrix. The transformation of the variables 7 a/ 3 defined by 
the matrix above is the 7 Q/ 3 = M 11 a M v ^ ^ . This transformation reads analytically: 



7n 

7l2 
713 

722 

723 

I 733 



e 2si 7n + s 4 (2e Sl 7 12 + s 4 72 2 ) 

e Sl s 2 7ll + (e Sl+S5 + s 2 s 4 ) 712 + e S5 s 4 722 

e Sl s 3 7n + (e Sl s 6 + s 3 s 4 ) 712 + e Sl 713 + s 4 s 6 722 + s 4 723 

e 2s5 722 + 2e S5 s 2 7i2 + S2 7n 



(5.7) 



,*r, 



723 + e S5 s 6 722 + s 2 713 + (e S5 s 3 + s 2 s 6 ) 712 + s 2 s 3 711 



733 + 2 s 6 723 + Sg 722 + 2 s 3 713 + 2 s 3 s 6 712 + sj j u 



The aforementioned equations define a transformation group G r of dimension r = 
dim(Aut(IV)) =6. In exact analogy to the previous Type II case (see (14.7J) ) the six 
group generators are found to be: 



X 1= 2 7ll 



d 



712- 



5 



+ 7i3- 



d 



0711 ^712 <^7i 3 

d d d 

X 2 =711-5 ^ 2 7i25 ^713-5 — 

0712 ^722 »723 

9 9 „ d 

X 3 =7u « ^7125 K2 7i3- 



^713 

X 4 =2 7l2 



'^723 

+ 722" 



+ 723" 



'^733 

5 



(5.8) 



X 5 =712- 



0711 ^7i2 c/713 

d d d 

+ 2 722 5 h723" 



' #712 ^722 ^723 

d d d 

X% =712-5 V 7225 V 2 723- — 

0713 «723 (^733 

The algebra g r corresponding to the group G r has the following non zero commutators: 



[X 1 ,X 2 }=X 2 , [X U X 3 ] 

[X 2 , X4] — X\ — X 5 , [X 2 , x 5 \ 
[^3,^4] = —Xe, [X4, X 5 ] 



x 3 , 


[Xi,X 4 ] 


— —Xi, 


x 2 , 


[X 2 ,X 6 ] 


= x 3 , 


— X4, 


[^5,^6J 


= x 6 



(5.9) 
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It is evident from the aforesaid commutators ( I5.9J) that the group is non abelian, 
hence we cannot diagonalize all generators at the same time. The derived algebra of g r 
is 

g r , = {[X A , X b ] : X A , X B E g r } => g r , = {X x - X 5 , X 3 , X 6 , } (5.10) 

the second derived algebra is 



g r „ = {[X A ,X B ] : X A ,X B Gg r ,}4 g r „ = {X 3 ,X 6 } 
and, finally, the third derived algebra is 

g r „, = {[X A ,X B ] : X A ,X B G g r „} =► g r „, = {0} 



(5.11) 



(5.12) 



Thus, the group G r is solvable since g r „, is zero. 

The relations (J5.9P imply that X 3 , X 6 commute with one another, as well as with the 
generator Y 2 , so that we can bring them to their normal form at the same time. To this 
purpose, the transformation of 7 a/ j is: 



7ii = e Ul 
712 = e Ul u 2 



7i3 = e" 1 (u 2 «3 + u s ) 



722 = e" 1 M 4 



(5.13) 



723 = e Ul (U3U4 + U2M5) 
k 733 = e Ul (-2 e m +u 2 3 u 4 + 2u 2 u 3 u 5 + uf) 

In these coordinates, the generators Y 2 ,Xj take the form 



Y 2 



d 
dui 



X 3 

d 



_d_ 



Xa 



X 4 = 2u 2 — +(-2u 2 2 + u 4 ) 

au\ ' ou 2 



d 

du 3 

d d d 

«5 — 2u 2 u 4 



2u 2 



X 5 = u 2 
X x = 2 



du 3 ' du 4 

d d d d d d 

a U3 T> V^u A —— X 2 = h 2 u 2 u 3 -^— 

OU 2 OU 3 OU4 OU 2 OU4 ou 5 

d d d d d 

u 2 — 2 m 4 7^ u 5 — 2 



d 

du 6 



(5.14) 



du\ du 2 dui du<a due 

i.e., the functions u%, u 3 , u§ appear in the field equations only through their derivatives. 

► The solution space in the gauge N a = 

Before we begin with the solution of the Einstein equations, we must examine which 
values are permissible for the functions Uj, i — 1, . . . , 6. The determinant of the matrix 

la/3 is: 

det[j a/3 ] =2e 3ui+U(i (u 2 2 - w 4 ) (5.15) 

so, the inequality u\ — U4 > must hold. 
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The two linear constraints E\ = 0, E 2 = in the variables (I5.13P give: 

E 1 =0 =>■ u 2 u 3 + u 5 = 

£?2 =0 =>■ «4 W3 + U 2 U5 = (5.16) 

The above system for u$, u§ has non zero solutions only in the non permissible case 
u 2 — W4 = 0, so that we obtain 

u-i = fa, u 5 = fa (5.17) 

However, these values of u 3 , u§ render 713,723 functionals depending on 711,712,722, see 
( 15.1 3p . Thus, it is possible to set 713, 723 to zero by a proper automorphism matrix (15. 6p . 
so that we have the following theorem: 

The general solution of Bianchi Type V corresponds to a block diagonal form in the 
basis of the 1-forms cr a . 

Having proved that the matrix 7 Q/ g is block diagonal with no loss of generality, we 
can repeat the former procedure for determining the generators, i.e., starting from the 
matrix: 

(7n 7i2 \ 
712 722 (5.18) 

733/ 

with the corresponding automorphism matrix 

/e Sl s 2 
M = s 4 e S3 I (5.19) 




In this case, the three generators Y 2 ,Xi + X 3 ,X 2 commute so that they can be diago- 
nalized at the same time. By applying the transformation: 



gUl+2«6 — 5 "4 


pUi+2v,6— 2 U4 ., 





„U1 +2 U(j — 2 U4 q. 


e «i+2« 6 -i U4 ( e «4 + u 2^ 











e Ul 



<y a p= \e Ul+2uG -2 U4 u 2 e Ul+2u «-2 UA (e U4 + ul) | (5.20) 

the generators take the form 

Y 2 = ^-,X 2 = ^-,X 1 + X 3 = ^- (5.21) 

OUi OU 2 OUq 

which means that the field equations are a first order system for the variables ill, ii 2 and 
iiQ. Consequently, the only variable the second derivative of which appears in the field 
equations is w 4 . In order to further simplify this situation, we will also use the generator 

v - d 
y i - at- 
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The variable u^ that does not appear in the generators Y 2 , X 1 , X 2 , can be used to 
parameterize time t, i.e., to apply the transformation 

t ->■ u 4 (t) = s, ui(t) ->■ Wi(s), w 2 (*) -> «2(s), %(*) ->■ «e(s) (5.22) 

Of course, this transformation can be applied only when the function u^{t) is not 
constant. So, we have to examine separately two cases, i.e., u^(t) being constant and 
not constant. 

5.1 The case u±(t) = const. 

Assuming that u 4 (i) = k&, the linear constraints E\ = 0, E 2 = are identically satisfied. 
The third linear constraint is transformed as 

E 3 = => u 6 = =► u 6 = k 6 (5.23) 

Substituting in the quadratic constraint ^0 = yields the following lapse function: 

N 2 = — e~ ki+ui (3 e ki ii\ - u 2 2 ) (5.24) 

Using this lapse function in the equations of motion f l2.5|) . we have as unknown functions 
Mi, u 2 . Returning to En = 0, it is evident that it is proportional to the quantity 

ii 2 (6 e kA u\ ii 2 - 2 ii\ + 3 e kA ii 2 ui-3 e ki ii x u 2 ) 

The vanishing of the quantity in the parentheses leads to a zero lapse, therefore it is 
rejected. The case u 2 = k 2 satisfies all field equations and gives for the lapse function 

N 2 = ^e Ul iil^(Ndt) 2 = (de u ^ 2 ) 2 => (N dt) 2 = (dr) 2 (5.25) 

Thus, choosing the time gauge U\ = 2 In r and the automorphism matrix 

( e fc 4 /4-fc 6 _ e -fc 4 /4-fc 6 U q\ 
e -fc4/4-fc 6 o (5.26) 

l) 

we obtain the metric 

d s 2 = _ dT 2 + T 2 (o .l )2 + r 2 (o . 2) 2 + r 2 (<T 3 )2 (5 _ 27) 

describing a flat spacetime. 
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5.2 The case u±(t) ^ const. 

In this case, putting u±(t) = hit 2 implies that the first two linear constraints are identi- 
cally zero, while the third linear constraint £'3 = yields 

E 3 = =► m 6 = =► u G = k 6 (5.28) 

The quadratic constraint E = defines the following lapse function: 

The substitution of the aforementioned lapse function A^ 2 and Uq in equation En = 
has as a result that the coefficient of u 2 is proportional to the quantity 

ii 2 (tiii - 1) . 

Consequently, before we solve equation En = for u 2 , we must secure that the above 
quantity is not zero. The vanishing of the term t ii\ — 1 leads to metric with signature 

( h+), hence it is non permissible. However, the vanishing of the term u 2 leads to a 

solution of the field equations as follows: Substituting u 2 = k 2 in equation E n = 0, the 
latter becomes proportional to 

(3tiii- l)(3t 2 Mi-3t 2 w 2 + 3tui + l) (5.30) 

The vanishing of the first term leads to incompatibility, while the vanishing of the second 
term gives 

u\ = fci2 + In I csch I kn H — j=- hit J J . (5.31) 

As a result, all field equations are satisfied. Now, using the transformation t = e^ 3( - T ~ kll \ 
redefining the constant k\ 2 = Inn 2 , and acting with the automorphism matrix 



M 



/ v3 h\\ , v3kii I,* \ 

e 2 B e — e 2 K& k 2 0\ 

v3fcii , 

e^r 1 "^ 



(5.32) 



V 01/ 

we obtain the metric 
d s 2 = k 2 ( - csch 3 r d t 2 + e~^ T cschr(<r 1 ) 2 + e^ T cschr(o- 2 ) 2 + cschr(cr 3 ) 2 ') (5.33) 

This metric does not admit a homothety, hence the constant k cannot be absorbed. The 
metric was first derived by Joseph [22J. 

Since f)5.30p does not vanish, we solve equation En = for u 2 and, by substituting 
in £"22 = the latter becomes proportional to the quantity 

(tu 2 + (3tui- l)u 2 ) (3t 2 Mi-3t(w 2 - 1)m 1 + m 2 -3£ 2 m 2 + l) (5.34) 
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The vanishing of the first term 

u 2 = e - 3ui k 2 t (5.35) 

leads to a solution of the field equations as follows: Solving En = for u\ and substi- 
tuting in Eu = we get k 2 = 0. As a result, equation Exx = yields 



• 2 1 • 1 

ux = u 1 --u x - — ^u x 



k 12 + \\x (sech I k u + — kit J j (5.36) 



that satisfies also the rest of the field equations. Now, using the transformation t h- >■ 
e -V3(5+fcn)^ Xl _^l (y3 ^ii — x), redefining the constant fc 12 = In A 2 , and acting with the 
automorphism matrix 

/ e~ k6 0\ 
M = e" fc6 (5.37) 

V l) 

we obtain the metric 

d s 2 = A 2 Q sech 3 £ d£ 2 + e-^ H sech£ (o- X ) 2 + e^ sech£ (° 2 ? + sech£ (<r 3 ) 2 ) (5.38) 

that does not admit a homothety, hence the constant A cannot be absorbed. This metric 
has a Euclidean signature. 

After having examined the case of the vanishing of the first term in f)5.34p . we proceed 
to the vanishing of the second term and get the following system of differential equations: 

ii 2 = - (u\ + 1) «2 (5.39a) 

v 

u x = — (3t 2 ul-3tui + 3tu x ul-ul- 1) (5.39b) 

3 1 

Equation f)5.39al) can be readily integrated since it is of first order for ii2{t) or, more 
easily, if we multiply by u 2 and put u 2 (t) = U(t), i.e. 

U2 u 2 = - (ii 2 2 + 1) ii 2 2 => - U = - (U + 1) U => 

dU 2dt , U 

In— =2 lnt + 2A;2i => 



U(U + 1) t U+l 

e 2fc 2 i+2 

U = — 2 =* u 2 = k 22 ± e- k21 V / l-e 2fc 2it2 (5.40) 

Substituting in (15.39bl) we have for U\. 

Ul = k 12 - In ( sinh (kn + —j= arccoth-\/l - e 2k2lt ' 2 J J (5.41) 

With these results, we unexpectedly obtain again Joseph's metric f )5.33p . 
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6 Bianchi Type IV 

For this type the structure constants are 

ri\ ri\ f~i2 f~i2 1 

u 13~ u 31 ~ u 23 ~ u 32~ 2 

C\ 3 = ~C\ 2 = i (6.1) 

C a p = for all other values of a[3^ 

Using these values in the defining relation (12. 2p of the 1-forms erf we obtain 

e~ x xe~ 
erf- ( -e" x | (6.2) 

The corresponding vector fields ££, (satisfying [£ Q ,£,g] = CP a a^) with respect to which 
the Lie Derivative of the above 1-forms is zero are: 

ii = -d y 6 = 9* Z,3 = d x + (y - z)d y + zd z (6.3) 

The Time Depended A.I.D.'s are described by 

(P{t) P{t)HcP{t)) x(t)\ 
A<* = P(t) y(t) (6.4) 

\ 1 / 

and 

x(t) (l*||y) -!/(*), y(t) (^f^y)''-^^))') ( 6 - 5 ) 

where P(t),x(t) and y(t) are arbitrary functions of time. As we have already remarked 
the three arbitrary functions appear in P a and thus can be used to set the shift vector 
to zero. 

The remaining symmetry of the EFE's is, consequently, described by the constant 
matrix: 

/e Sl s 2 s 3 \ 
M = e Sl s 4 (6.6) 

\0 l) 

where the parametrization has been chosen so that the matrix becomes identity for the 
zero value of all parameters. 

Thus the induced transformation on the scale factor matrix is 7 a/ 3 = M£Ma / -f flu , 



26 



which explicitly reads: 

7n = e 2si 7n 

712 = e Sl s 2 7i2 + e 2si 7i2 

7i3 = e Sl (s 3 711 + s 4 712 + 713) 

722 = e 2 Sl 722 + 2 e si s 2 712 + s 2 7n 

723 = e Sl {s 3 712 + s 4 722 + 723) + s 2 (« 3 7n + 7i3 + s 4 712) 

k 733 = -5 3 2 711 + S3 (712 + 713) + S4 2 722 + 2 S 4 (S3 712 + 723) + 733 

The previous equations, define a group of transformations G r of dimension r 
dim(Aut(IV)) = 4. The four generators of the group are: 



(6.7) 



9 n d d n d d 

X ! = 2 7u- h 2 7i2^ h 7i3^ h 2 722 h 723 — 

»7ll O712 C?7i3 722 723 



x 2 

X 3 = 711 



d n d d 

7iio ^ 2 ^ ^7i3 



'#7l2 

d 



Xa 



712- 



<97l 3 

<9 



+ 712 
+ 722 



^722 

d 

C?723 
9 



+ 2 713 
+ 2723 



^723 

d 



^733 

<9 



(6.8) 

(6.9) 
(6.10) 
(6.11) 



»7l3 <?723 <^733 

The algebra g r that corresponds to the group G r has the following table of non-zero 

commutators: 

[Xi,X 2 ] = 0, [X l5 X 3 ] = X 3 , [Xi,Xa] = X 4 , , . 

[X 2 ,X 3 ]=0, [X 2 ,X 4 ]=X 3 , [x 3 ,x 4 ] = o l °' iZJ 

As it is evident from the above commutators (I6.12p the group is non-abelian, so we 
cannot diagonalize at the same time all the generators. However, if we calculate the 
derived algebra of g r , we have 



g r , = {[X A ,X B ] : X A ,X B E g r } =► g r , = {X 3 ,X 4 } 
and furthermore, it's second derived algebra reads: 

g r „ = {[X A ,X B ] : X A ,X B Eg r ,}^ g r „ = {0} 



(6.13) 



(6.14) 



Thus, the group G r is solvable since the g r „ is zero. As it is evident X 3 ,Xa,Y 2 
generate an Abelian subgroup, and we can, therefore, bring them to their normal form 
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simultaneously. The appropriate transformation of the dependent variables is: 



7ii = e ni 



712 = e Ul u 2 



7i3 = e Ul (-e" 4 u\ + u 3 + u 2 u 5 ) 



(6.15) 



722 = e ni+n4 



723 = e ni (w 2 (w 3 - 1) + e" 4 w 5 ) 

k 733 = e Ml (e- 2M6 + w 3 2 - e~ U4 w 2 (2 u 3 - 1) + 2 w 2 u 5 (u 3 - 1) + e M4 w 2 ) 
In these coordinates the generators Y 2 ,Xa assume the form: 

~Y d v d v d 

2 dui ■* du$ 3 q Us 

X 2 = ^ + (e-- M 2-«5)^ + 2e-- M 2^ + e- 2 -(e--2 M 2 )^ (6.16) 

*i = 2 A-«3^ + (e-« 2 -«5)^ + i 

Evidently, a first look at (I6.15P gives the feeling that it would be hopeless even to 
write down the Einstein equation. However, the simple form of the first three of the 
generators (16.161) ensures us that these equations will be of first order in the functions 
Mi, ii 3 and U5. 

6.1 Description of the Solution Space 

Before we begin solving the Einstein equations, a few comments for the possible values 
of the functions Ui, i = 1, . . . , 6 will prove very useful. 
The determinant of j a /3, is 

det[ laP ] = e 3 " 1 " 2 " 6 (e U4 - M 2 ) (6.17) 

so we must have e Ui > u\ . 

The two linear constraint equations, written in the new variables (16.151) . give 

E 1 = => e Ui u 3 + u 2 2 u 4 + u 2 (e U4 u 5 - u 2 ) = (6.18) 

E 2 = => (3 e" 4 + u 2 ) ii 2 - e Ui (1 + 3 u 2 ) ii 3 - (u 2 + 3 e" 4 ) (u 2 w 4 + e" 4 u 5 ) = (6.19) 
Solving this system for the functions u 3 ,ii5 we have 

u 3 = 0, ii 5 = e~ Ui (u 2 - u 2 u 4 ) (6.20) 
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yielding to 

u 3 = k 3 , u 5 = k 5 + e"" 4 u 2 (6-21) 

Now, these values of 113,115 make 713,723 functionally dependent upon 711,712,722 (see 
(I6.15P ). It is thus possible to set these two components to zero by means of an appropriate 
constant automorphism. 

Without loss of generality, we can start our investigation of the solution space for Type 
IV vacuum Bianchi Cosmology from a block- diagonal form of the scale-factor matrix 
(and, of course, zero shift) 

flu 712 \ 
7«/3 = 712 722 (6.22) 

V 733/ 

Note that this conclusion could have not been reached off mass-shell, due to the fact that 
the time-dependent Automorphism (16 .4p does not contain the necessary two arbitrary 
functions of time in the (13) and (23) components (besides the fact that all the freedom 
in arbitrary functions of time has been used to set the shift to zero). As we have earlier 
remarked, since the algebra (16.12p is solvable, the remaining (reduced) generators X\, X 2 
(corresponding to diagonal constant automorphisms) as well as Y 2 continue to define a 
Lie-Point symmetry of the reduced EFE's and can thus be used for further integration 
of this system of equations. 

The remaining (reduced) automorphism generators are 

^ d d d 

X\ = 711 7; h 712 7; h 722 



#711 ^7i2 ^22 

X 2 = lii 7 ^— + 2 ll2 —- (6.23) 

0712 0722 

The appropriate change of dependent variables which brings these generators -along with 
Y 2 - into normal form, is described by the following scale-factor matrix : 

(g«l+2u 6 g«i+2u6 u Q \ 

e Ul+2ua u 2 e Ul+2ue {u\ + u 4 ) (6.24) 

e Ul J 

The generators are now reduced to 

d d d 

Y 2 = —,X 2 = —,X 1 = — (6.25) 

au\ ou 2 ouq 

indicating that the system will be of first order in the derivatives of these variables. 
The remaining variable u 4 will enter, (along with it 4 , H4) explicitly in the system and is 
therefore advisable (if not mandatory) to be used as the time parameter, i.e. to effect 
the change of time coordinate 

t -)■ u 4 (t) = s, u x {t) ->■ ui(t(s)), u 2 (t) ->■ u 2 (t(s)), u 6 (t) -> u 6 (t(s)). (6.26) 

This choice of time will of course be valid only if u 4 is not a constant. We are thus led 
to consider two cases according to the constancy or non-constancy of this variable. 
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6.1.1 Case I: u 4 (£) = fc 4 

In these variables the first two linear constraint equations are identically satisfied, and 
the determinant of j a p, is 

detfrap] = e 3ui+4u6 h 

so we must have fc 4 > . The third linear constraint reads 



E 3 = => u 2 + 4 k 4 ii 6 = => u 2 = k 2 - 4 k 4 u 6 (6.27) 

Substituting this value into the quadratic constraint equation E we obtain for the lapse 
function 

N2 = 10 / fc4 , 1 eUl (3 ^? + 8 «! ^ 6 - 4 (4 k 4 ~ 1) ul) (6-28) 

iZ fc 4 + 1 

If we substitute this value of the lapse function into the equations of motion (12.51) . we 
are left with the unknown functions Ui,Uq. The strategy we follow is to solve one of 
(12. 5p for a second derivative of some function, say Uq and replace the result into the rest 
of the equations. In order to do that, we must ensure that the coefficient of Uq does not 
equals to zero. Looking at En = the coefficient of Uq is proportional to 

iii (ill + 2 (4 fc 4 + 1) u 6 ) =>■ {ui — k\ or U\ = k\ — 2 u 6 — 8 fc 4 u 6 } 

Thus we are forced to examine the above equalities before we solve E u = 0, for uq. 

The first possibility U\ = k\ yields an inconsistency, so we are left with the second, 
i.e. 

Ml = k\ — 2 Uq — 8 fc 4 Uq (6.29) 

The above choice satisfies all the spatial equations E a p = and gives the lapse function 

Redefining the constants k\ = ln/t 2 , k 4 = 4( -~^ 1 - ) , choosing a time parametrization 
uq = (« — 1) ln(t), and using the automorphism matrix (16.61) with entries s\ = 0, s 2 = 
—k 2 , S3 = 0, s 4 = we arrive at the line element 



ds 2 = -« 2 (dt) 2 + t 2fl (a 1 ) 2 + 2£ 2 > Into- 1 



<x 2 



+ t 2 ^ Utft-j^-^] (<r 2 ) 2 + £ 2 (<r 3 ) 2 (6.31) 

with t > 0, < fi < 1. In the above line element we have dropped the constant k since 
this line element admits a homothetic vector field 

H = d t + (fxz-y(fx-l))d y ~ z(fx - 1) d z (6.32) 
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Line element (16.3 ip was first derived by Harvey and Tsoubelis [23] and admits, besides 
the three killing fields ( 16 .3p three more, namely 

(6.33) 
d x + fxd v + f 2 d z (6.34) 



£4 


=e 


-x/v df+ e 


-x/n 

t 


d x 


6 


=y 


e- x/ » dt + 


t 


x/fi 


c 




„-x/n Q , 


ze~ 


x/fj, 



d x + f 2 d y + f 3 d z (6.35) 



with 



fi= * \_ 2fi+1)3 ( 4 ^ ( 2 A*- I) 2 (Ai-1) In 2 1 

-8 y u(/i-l)(2/i- l)(-/i+(2/i-l)x) hit 

- /i + 4 (// - 1) (/i + (1 - 2 /i) x) 2 



/ 2 = ^ (-2^+1)2 (^(2^-l)lnt + / i+(l-2 / i)x) 

, 4//( / /-i) e (^- 1 )^r 2 ^ 1 

/3 = ^— T^ 6 - 36 

— 2//H- 1 

There is thus a G§ symmetry group acting (of course, multiply transitively) on each V 3 
of this metric. However, it is interesting to note that we have not imposed the extra 
symmetry from the beginning, but rather it emerged as a result of the investigation 
process. 

Having ensured that the coefficient of Uq at En = is not zero, we can solve this 
equation for u e and insert the result into the rest of the spatial equations. But doing that 
we end up with a zero lapse, indicating that the only solution for this case is described 
by the line element ( 16.3ip . 

6.1.2 Case II: u A {t) = t 

In this case the determinant of the sale factor matrix equals to 

det[ la p\=e^ +iu H 

so we must demand that t > in order for ^ a p to be positive defined. 

The first two linear constraints are identically zero while the third one E 3 = can 
be used to define the function u 2 



-t-lj; 



E 3 = 0^u 2 + 4:tu 6 + l = 0^U2 = k2-t-4: tu 6 dt (6.37) 
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and the quadratic constraint E = defines the lapse function N 2 



j u 1 



N 2 



12* + 



- (3tul + 8tuiU 6 -U(U- 1) ii 2 6 - 2 (At - 1) ii 6 + 2u ± - l) ( 6 - 38 ) 



Substituting the above values of the lapse N 2 and the function u 2 in equation E 33 = 
we find the coefficient of Uq is proportional to 

Mi (Atu x -{At- 1) (Atu 6 + 1)) 

a quantity that can be safely regarded different from zero, since it's nihilism leads either 
to zero lapse or to inconsistency. Thus we can solve E 33 = for uq and substitute it to 
En = 0. In order to solve this equation for v,i we must be assured that it's coefficient 
does not equals to zero. Setting this coefficient equal to zero we arrive to the following 
equation 

«i = l + (4t- l)u 6 

which is unacceptable because it leads to inconsistency. After solving equation En = 
for iii we finally arrive to the following polynomial system of first order in iix, w 6 



iii = (ui\Bi\u 6 ), uq = (ui\B 2 \u 6 ) 



(6.39) 



where we have used the notation (iii\ 
matrices Bi, B 2 given by 



luiu 2 uf) and \iiij = (nj| 4 with the 4x4 



Bi 



i(12t+l) 



/— 4t — 1 -32t 3 + 2t -64t 3 + 4t 0\ 

24t 2 + l lQt 2 + At 8t 2 (4t+l)(12t- 1) 

12t 2 -t -16 1 2 

\ -6t 2 0/ 



B, 



i 



2t(12t+l) 



/12 1 + 3 144£ 2 + 8£-6 16 1 (36 1 2 + It - 1) 16 1 2 (48 1 2 + 8t - 1)\ 

-4 -24 1 -32 1 2 

-6t -12t 2 

\ / 



Due to the form of Bx, B 2 (their components are rational functions of the time t), system 
(I6.39P can be partially integrated with the help of the following Lie-Backlund transfor- 
mation 



u-i 



u 6 



(48 1 2 + 16 1 + 1) r(t) - 2 (12 1 - 1) tan r(t) 
4v^(12t+l) 

\/i(12t+ 1) r(t) + QVt tanr(t) + 24 1 + 2 
8£(12£+1) 



(6.40a) 
(6.40b) 
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yielding the single second order ODE for the function r(t) 

(6t+l) tanr + Gv^ . 



ftanr + Vt)r 2 -2 ^ t+ ^ tanr + 
V / V^(12i+1 



36 1 2 tan 2 r + 36v / t3tanr-12t-l . 

+ ■=- - (6.41) 

V¥(12t + 1) 2 

This equation contains all the information concerning the unknown part of the solution 
space of the Type IV vacuum Cosmology Unfortunately, it does not posses any Lie-point 
symmetries that can be used to reduce its order and ultimately solve it. However, its 
form can be substantially simplified through the use of new dependent and independent 
variable (p,u(p)) according to r(t) = =|= arcsin y= , t = 6( - 1 _ 1) thereby obtaining the 

equation 

l-« 2 -2 (1-M 2 ) 2 , nAn , 

u = + - — => u 2 = — \ , - ; (6.42) 

y/6(p-l)(f?-u*-l) 6(p-l)(p 2 - M 2 -l) l ; 

This equation is a special case of the general equation 

2 (1-U 2 ) 2 , > 

U = (. + Xp)(p 2 -u 2 -l) (6 - 43) 

with the values k = —6, A = 6. The general solution of (16.431) was first given in [24] and 
can be obtained as follows: We first apply the contact transformation: 



8 _. 4(2f-l) .... k 4 ,,.. 

u(p) = -j 2/(0 + -^ — - y'(Z) p = -j + x y ^ 

u(p) = 2£-l u(p) 



(6.44) 



2y"(0 

which reduces it to 

e(Z-i) 2 y" 2 = -4y'(Zy'-y) 2 + 4y' 2 (Zy'-y)-^y' 2 + ^^y' (6.45) 

This equation is a special form of the equation SD-Ia, appearing in [25], where a 
classification of all second order second degree ordinary differential equations was per- 
formed. The general solution of (I6.45P is obtained with the help of the sixth Painleve 
transcendent w := Pyi(a!, /3, 7, 8) and reads: 

e (j - 1) 2 ( , w(w-i) 

y = -a — 7 — 7T7 — 7^ w - 



iw (w - i)(w - 1) V. £(£-i 

1/, , .A7-X2/, n ..., P (, 2^ 



+ -(l ± V2 a ) 2 (l-2.)-^l w 
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where the sixth Painleve transcendent w := Pvi(cn, /3, 7, S) is defined by the ODE: 



w 



1/1 1 I 

+ - + 



2 \w — 1 w w — £ 
u> (w — 1) (w — £) 



/2 

W — 



1 1 1 

+ T + 



e-i e ™-e 



w 



e a - 1) 



a + f3 — + 7 



G-i) , ^(e-i; 



w 



., +5 



(6.47) 



The values of the parameters (a, (3, 7, 5) of the Painleve transcendent, can be obtained 
from the solution of the following system: 

a-/3 + 7-<5±v / 2^+l = ~ (6.48a) 



(/3 + 7) (a + 5±V2a) = (6.48b) 



k 2 -A : 



(- - »') (n - S± \/2d - I ) + - (^a-/3-7 + 5±v / 2^) ]( . 



I (7-/3) (a + 5±V2^) 2 H-^ (/3 + 7) 2 (a-<J±v^a + l) 







(6.48c) 
(6.48d) 



Inserting in f )6.48[) the values of « = —6, A = 6 for Type IV, we have twenty-four 
solutions (counting multiplicities) of this system. In order for the parameters (a, ft, 7, 5) 
to be real numbers we end up only with three possibilities 



1 3 3 1 
2'~2'2' 2 



01,13,7,6) = 

a, p,y,8) = (2 + ^,0,0,-1 
a,/3, 7 ,5) = (2-^,0,0,-1 



(6.49a) 

(6.49b) 
(6.49c) 



Gathering all the pieces the final form of the general line element describing Bianchi 
Type IV vacuum Cosmology is 



ds 2 



3«l(£) 



4£(e-i; 



(dO 2 + v / OPW(^ 1 ) 2 + 2 v^-iMO^O o- 1 



+W^p ( 4 ^(0 |/(0 + 1) (- 2 ) 2 + e" 1 ® (- 3 ) 2 



(6.50) 



where 



l-2e-2y(Q 



y(0 



2e(e-i) 



2^(0(^-1) 



(6.51) 



and y(£) is given by (16.461) . 
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7 Discussion 

The present work completes the first phase of the programm initiated in [10], utilizing the 
automorphisms of the various Bianchi Cosmological Models as Lie- Point Symmetries of 
the corresponding Einstein's Field Equations with the aim of uncovering their solution 
space. The power of the method lies in the fact that it constitutes a semi-algorithm 
which, if successfully applied, results in the acquisition of the entire space of solutions. 
This successful implementation had, so far, been carried out in the case of Bianchi Type 
III ([21]) and Vllh ([26]; while the present communication covers Types /, II, IV and 
V. In all cases considered where the general solution is expected to have three essential 
constants (III, IV and Vllh), it is given in terms of the sixth Painleve transcendent 
Pvi, along the way with all the known particular solutions. It is noteworthy that these 
known metrics are rediscovered in a systematic way and without any extra assumption, 
in contrast to how they were originally obtained. The case of Types J, II, and V 
is characteristic: The general solutions, not considered as such at the time of their 
first derivation and containing 1 or 2 essential constants, were produced with the aid 
of various simplifying ansatzen in a time scale of half a century; Kasner (1921), Taub 
(1951), Harrison (1959), Joseph (1966). Here, they are comprehensively re-acquired along 
with the solutions not attributed to any one else which, to the best of our knowledge, 
are new: The Type I metric (13. 32ft . the Type II metrics (14.281) (Euclidean signature), 
(I4.3?p . (I4.5ip (Euclidean signature) and the Type V metrics (I5.27P and (I5.38P (Euclidean 
signature). The production of metrics with Euclidean signature may, at first sight, strike 
as odd; since our staring point is a line element of Lorenzian signature. However, it is 
made possible by allowing the lapse to be determined through the quadratic constrained 
equation instead of prescribing it by an ab initio choice of time gauge. 

In the remaining Bianchi Types VIII, IX and the exceptional Vlh the number of 
existing automorphisms is not sufficient to allow our method to reduce the problem to 
a single second order ODE in one unknown function, but rather to a third order one. 
We strongly suspect it to be an equivalent form of the Chazy type equations. However 
the task of proving it involves the search for an appropriate Lie-Backlund transformation 
which is highly non-trivial and non-algorithmic. We plan to return if and when there is 
something concrete to be reported. 

Some directions for future work include the application of the method in the presence 
of matter sources and/or in higher dimensions. 
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